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Some problems related to construction of the e-expansion of dimensionally regulated Feynman integrals are 
discussed. For certain classes of diagrams, an arbitrary term of the e-expansion can be expressed in terms of 
log-sine integrals related to the polylogarithms. It is shown how the analytic continuation of these functions can 
be constructed in terms of the generalized Nielsen polylogarithms. 
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1. Dimensional regularization jl| is one of the 
most powerful tools used in loop calculations. In 
some cases, one can derive results valid for an ar- 
bitrary space-time dimension n, usually in terms 
of various hypergeometric functions. However, 
for practical purposes the coefficients of the ex- 
pansion in e are important, where the regulator e 
corresponds to the difference between n and the 
(integer) number of dimensions of interest. Below 
we shall usually imply that n — 4 — 2s. In multi- 
loop calculations higher terms of the e-expansion 
of one- and two-loop functions are needed, since 
they may get multiplied by poles in e, not only 
due to factorizable loops, but also as a result 
of applying the well-known reduction techniques 

In refs. ftUq], it was shown that the log-sine 
integral functions (see, e.g., in p|, chapter 7.9), 
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happen to be very useful to represent results for 
higher terms of the e-expansion. 

For instance, for the one-loop two-point func- 
tion j( 2 )(n; vi, V2) with external momentum k, 
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masses mi and m.2 and unit powers of propaga- 
tors, the following result for an arbitrary term of 
the e-expansion has been obtained in pjfy : 
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where 

cosT^ = (to 2 - to 2 + fc 2 )/(2mi\/fc2) 



(2) 



COST 02 
COS T12 



(ml - m\ + fc 2 )/(2TO 2 vfc2) 

(TO 2 + TO 2 — fc 2 )/(2TOlTO2) , 
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whereas the "triangle" function A is defined as 

A(x, y, z) = 2xy + 2yz + 2zx - x 2 - y 2 - z 2 . (4) 

One can see that 7"i2 + '7"oi+' r 02 = ^- In fact, these 
angles can be associated with a triangle whose 
sides are mi, TO2 and Vfc 2 - Moreover, the area of 
this triangle is |a/A(to 2 , to 2 ,, fc 2 ,). For details of 
geometrical description, see in J?J. 

Note that the values of Lsj(7r) can be ex- 
pressed in terms of Riemann's C function, see 
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Eqs. (7.112)-(7.113) of §. The infinite sum with 
Lsj(7r) in (|3j) can be converted into T functions, 
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by the imaginary (and real) parts of the general- 
ized Nielsen polylogarithms (see, e.g., in poj ), 

(_ 1 \a+b-i r \ n a - 1 n n b (l-z£) 



2. The e-expansion (0) is directly applicable 
in the region where A(mj, m 2 , k 2 ) > 0, i.e. when 
(mi — ra 2 ) 2 < k 2 < (mi+m 2 ) 2 . In other regions, 
the proper analytic continuation of the occurring 
Lsj (9) should be constructed. To do this, it is 
convenient to introduce the variable (cf. in ref. I| ) 



ln(— z) = ln(z) — i(77r, 
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where the choice of the sign a = ±1 is related to 
the causal "+i0" prescription for the propagators. 
Since Lsi(0) = —0, we get 



ia[Lsi (tt)-Lsi (9)} = ln(-z). 



(7) 



For the next order, we can use the fact that 
Ls 2 (0) = Cl 2 (0), where (see in §) 



01,(9) = 



-Li,- 



j even 



i+Lij (e- i9 )] , j odd 



(8) 



is the Clausen function, whereas Lij is the poly- 
logarithm. In other words, C\j(9) corresponds 
either to the imaginary part or to the real part of 
Li, (e l6 ) , depending on whether j is even or odd. 
Therefore, the analytic continuation reads 

i<7 [Ls 2 (vr)-Ls 2 (0)} = -\ [Li 2 (z)-U 2 (1/z)} , (9) 

where Ls 2 (tt) = 0. The result for the e-term of 
the two-point function was obtained in Q . 

To proceed further, we need similar relations 
between higher L,Sj (9) and the imaginary (or real) 
parts of the polylogarithms. For j = 3, Ls 3 (0) 
can be expressed in terms of the imaginary part 
of Li3 (l — e 10 ) . Then, the imaginary part of 
L14 (l — e l9 ) is already a mixture of Ls4(#) and 
CU(9), whereas its real part involves the general- 
ized log-sine integral Ls^ 1 ^). All these relations 
can be found in |6j (some misprints are mentioned 
in ID). However, attempts to generalize these re- 
sults to higher functions show that the relations 
get more and more cumbersome. 

Instead of going that way, we suggest to con- 
sider how the higher Lsj functions are generated 



where S a ,i{z) = Li a+ i (z) . We obtain 

Re S h2 (e ie ) = ±C1 3 (9) + ±Cs - \ (tt-0) Ls 2 (9) , 

Im S h2 (e w ) = -|Ls 3 (9) - ±9 (9 2 -3it9+3tt 2 ) , 
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+^ A + h^e-^ 2 9 2 ' 1 



n9 i 

48 192 



1 ai 



90 1 48 

Im Si, 3 (e i9 ) = iLs 4 (9) + ±C1 4 (9) - \ir(a 
+ l(n-9)C\ 3 (9)-±(n-9 fLs 2 (9), 

where the generalized log-sine integral (see in || ) 
is defined as 



Ls\ k) (9) 
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In particular, Ls^ " 1 (9) = hsj (9). Using these re- 
lations, we can express the Ls 3 and LS4 functions, 

iff [Ls 3 (tt) - Ls 3 (9)} = 

S 1<2 (z)-S 1>2 (l/z)-±hi 3 (-z), (12) 
iff [Ls 4 (tt) - Ls 4 (9)} = -3 [S h3 (z) - S 1<3 (l/z)} 

+ | [Li 4 (z) - Li 4 (1/z)] 
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-f [Lis (z)+Li 3 (l/z)] ln(-z) 
+|[Li 2 (z) - Li 2 (1/z)] ln 2 (-z) , 



(13) 



where Ls 3 (w) = —j^n 3 , LS4 (n) = §7r£ 3 . 

We have also constructed further expressions, 



ia [Ls 5 (vr) - Ls 5 (9)} = 12 [S 1A (z) - Si A (l/z)} 

-3[S 2 , 3 (z)-S 2>3 (l/z)} 

+3[Si t3 (z)+S lt3 {l/z)} ln(-z) 

+ ±ln 5 (-z), (14) 
\a [Ls 6 (tt) - Ls 6 (9)} = -60 [S h5 (z) - S h5 (l/z)} 

+15 [S 2A (z) - S 2A (l/z)} 

-f [Li 6 (z)-U 6 (l/z)} 

-15[S 1A (z) + S 1A (l/z)} ln(-z) 

+ f [Li 5 (z) + Li 5 (1/z)] ln(-z) 
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-f [Li 4 (z) - Li 4 (1/z)] ln 2 (-z) 
+| [Lis (z) + U 3 (l/z)} ln 3 (-z) 
-^[Li 2 (z)-U 2 (l/z)} ln 4 (-z) 



(15) 



with Ls 5 (tt) = -i|^ 5 , LsgfTr) = f Ks + f^Ca 
Substituting (0), ®, (HI)-® into Eq. @) 
(taking = Tq X or = Tq 2 , denoting the z's from 
Eq. (|^) as zi and z 2 , and setting cr = —1), we ar- 
rive at the analytic continuation of the terms of 
the e-expansion, up to order e 5 . In fact, we have 
also obtained results for higher Lsj functions, up 
to j = 10, which allowed us to reach the order e 9 . 

It is instructive to consider the limit to 2 — > 0. 
Introducing the variables x — m 2 /k 2 and y = 
m 2 /k 2 (and remembering that a = —1), we get 



y 
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In [A(m 2 ,TO 2 ,fc 2 )/(fc 2 ) 2 
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Then, S a j,(l/z\) can be converted into S a ,b(^i) 
by means of known relations given in Ref. pOj . 
After this, the limit y — » (7712 — > 0) can be 
taken, since all In y terms are cancelled. The 
obtained expressions can be simplified by trans- 
forming S , a ) b(z 2 ) into S a ^{l/z 2 ). In such a way, 
we also avoid appearance of terms like S a ,f>(—1). 
Note that for this limit the terms up to order e 3 
can be extracted from Eq. (A. 3) of ref. Our 
expressions are in agreement with their results. 
In fact, using hypergeometric representation 
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(see, e.g., Eq. (10) of p2| ), an arbitrary term 
of the e-expansion can be obtained. Employing 
Rummer's relations for contiguous functions, one 
can transform the 2F1 function from Eq. ([To) into 
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with u = k 2 /m 2 . This (transformed) 2-F1 func- 
tion can be expressed in terms of a simple one-fold 



parametric integral, 
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Expanding the integrand in e, we arrive at 
J (2) (4-2e;l,l) 
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i=i fc=i J 
which agrees with the results discussed earlier. 

3. In ref. [|| it was shown that similar explicit 
results can be constructed for the off-shell mass- 
less one-loop three-point function with external 
momenta pi, p 2 and p 3 (pi + p 2 + P3 = 0), 

J(t7; v 1 ,v 2 ,v 3 \pl,pl,p 2 .) 

dn 
T 

(18) 



[{p 2 -r)T [(Pi+r) 2 r(r 2 y^ 
as well as for the two-loop vacuum diagram with 
arbitrary masses mi, 7772 and 7773, 

I(n; v ll v 2l v 3 \rr? 1 ,m\,m 1 3 ) 

d p d <? 



(pa-m?)^ (g 2 -^ 2 ) 1 ' 2 [{p-q) 2 -m 2 ]^ ' (19) 

Results for general 77 and i>\ (in terms of hypergeo- 
metric functions of two variables) are available in 
Refs. JPJT^Q. According to the magic connec- 
tion |8|, these integrals are closely related to each 
other. For example, in the case v\ = v 2 = v 3 = 1 
this connection (see Eq. (16) of ||) yields 

J(4-2 £ ;l,l,l) = 7 r- 3 -i 1+2e {p 2 p 2 p 2 y e 

x r(i 1 -2 £ £ ) ) 7(2 + 2£;1,1 ' 1) ' (20) 

where we assume that pf <-*■ to 2 . Below we shall 
omit the arguments pf and to 2 in the integrals J 
and /, respectively. 

Then, using exact results in terms of 2-F1 func- 
tions Jl4] , p^[ [5[, in combination with the formula 

£fc^ L8j+l(2 „ = _ 2 „ r <i-^) 
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the following results have been obtained in || : 
J(4-2 £ ;l > l ) l) = 2 7r -M^ ni+Ws) 
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where the angles ^ (i = 1, 2, 3) arc defined via 

^i = (p^+pi-P?)/(2vSi), 



COS( 



COS l»2 



COS 



^2 = (P 2 +P 1 -P^)/(2 V /^ 2 ), 
03 = b 1 +^-p 2 )/(2vS 2 ) (24) 



(remember that p\ «-> mf for the integrals I), so 
that </>i + 02 + 03 = tt- Note that the angles 
9i from (lj^l are related to (pi as 0j = 2<^. By 
analogy with the two-point case (|J), the angles 
(pi can be understood as the angles of a triangle 
whose sides are , \Zp! an d x/pf, whereas its 
area is ±^A(p 2 ,p 2 ,p 2 ). 

For the two lowest orders (e and e 1 ), we re- 
produce eqs. (9)-(10) from [||. Useful represen- 
tations for the e° terms of both types of dia- 
grams can also be found in Jl(|. Moreover, in 
Eq. (26) of a one-fold integral representation 
for J(4— 2e; 1, 1, 1) is presented (for its generaliza- 
tion, see Eq. (7) of B). Expanding the integrand 



in e, we were able to confirm the e-expansion (|22| ) 
numerically. 

To construct the analytic continuation of the 
terms of the e-expansion (|2^) and ( 23]) , we need 
just to apply substitutions (Q), (|), (||)-(jl|), 
with 9 = 2<pi (i = 1,2,3). The remaining terms 
Lsj+2 (tt) can actually be treated in the same 
way, if we substitute 9 = 0. In any case, their 
values are known (in terms of £ function) and 
can be summed into a combination of T func- 
tions, Eq. (||). We get three variables Zi = e 2lCT( ^, 
see Eq. (^|), such that ziz 2 z 3 = 1. The causal 
prescription requires to take a = +1 for the J- 
integrals and a = —1 for the /-integrals. Using 
the substitutions presented above, we obtain the 
analytic continuation of the results ( p2| ) and ( |23| ) 
up to e 4 . We note that the result for the e-term 
was known jl7|,|| (in terms of LI3). 

When the masses are equal, mi = m 2 = ma 
m (this also applies to the symmetric case p\ = 



pi 



Pi 



p ), the three angles ■ 



7r/3, whereas A(m 



m , m 
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in this case the r.h.s. of Eq. (J 
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3m 4 . Therefore, 
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For instance, in the contribution of order e the 
transcendental constant Ls3(27r/3) appears. This 
constant was discussed in detail in gj. The fact 
that Ls3(27r/3) occurs in certain two-loop on- 
shell integrals and three-loop vacuum integrals 
has been noticed in []l8|-po||. Moreover, in Jjj| 
it was observed that the higher- j terms from ( |25| ) 
form a basis for certain on-shell integrals with a 
single mass parameter. Connection of Ls3(27r/3) 
with multiple binomial sums is discussed in |2lf| . 
We also note that in the constant Ls3(7r/2) 
appeared. 

4. One of the interesting problems is to con- 
struct terms of the e-expansion for the one-loop 
three-point function with general masses. In this 
sense, the geometrical description seems to be 
rather instructive. The geometrical approach to 
the three-point function is discussed in section V 
of (see also in Q]). This function can be rep- 
resented as an integral over a spherical (or hyper- 



5 



bolic) triangle, as shown in Fig. 6 of Q|, with a 
weight factor l/cos 1 " 2 ^ (see eqs. (3.38)-(3.39) 
of 0). This triangle 123 is split into three trian- 
gles 012, 023 and 031. Then, each of them is split 
into two rectangular triangles, according to Fig. 9 
of . We consider the contribution of one of the 
six resulting triangles, namely the left rectangu- 
lar triangle in Fig. 9. Its angle at the vertex 
is denoted as ^(pf 2 , whereas the height dropped 
from the vertex is denoted 7712. 

The remaining angular integration is (see 
eq. (5.16) of §) 



1 

2e 



tan 2 T] 12 s 
cos 2 (fi 



ff 2 /2 



2 ^(J + 1)! 

J-U 



dip hi 3 



+ 1 



tan 2 7712 
cos 2 ip 



.(26) 



First of all, we note that the l.h.s. of Eq. (26) 
yields a representation valid for an arbitrary e 
(i.e., in any dimension). To get the result for the 
general three-point function, we need to consider 
a sum of six such integrals. The resulting repre- 
sentation is closely related to the representation 
in terms of hypergeometric functions of two argu- 
ments [£3| (see also in |24j] for some special cases) . 

In the limit e — > we get a combination of 
CI2 functions, eq. (5.17) of 0. Collecting the re- 
sults for all six triangles, we get the result for the 
three-point function with arbitrary masses and 
external momenta, corresponding (at e = 0) to 
the analytic continuation of the well-known for- 
mula presented in j2j| . The higher terms of the e- 
expansion correspond to the angular integrals on 
the r.h.s. of Eq. (p6[). The problem of construct- 
ing closed representations for these terms, as well 
as their analytic continuation, is very important. 
We note that the e-term of the three-point func- 
tion with general masses has been calculated in 
|9] in terms of Li3 . 

5. We have shown that the compact structure 
of the coefficients of the e-expansion of the two- 
point function (J3), the massless off-shell three- 
point function (|22 ) and two-loop massive vacuum 
diagrams (p3|), in terms of log-sine integrals, al- 
lows to perform analytic continuation in terms of 



generalized Nielsen poly logarithms (10), in some 
cases ([[7]) even for an arbitrary order of the e- 
expansion. It is likely that a further generaliza- 
tion of these results is possible, e.g. for the three- 
point function with different masses, two-point 
integrals with two (and more) loops and three- 
loop vacuum integrals. In particular, numeri- 
cal analysis of the coefficients of the expansion 
of certain two-point on-shell integrals and three- 
loop vacuum integrals |l^| shows that in some 
cases the values of generalized log-sine integrals 
Ls^ , Eq. (|ll|), may be involved. For instance, 

in ref. [|l9| it was shown that Ls^ 1 "* (27r/3) is con- 
nected with V^i from pFjfl . 

The fact that the generalization of LS2 = CI2 
goes in the Lsj direction, rather than in Clj di- 
rection (see Eq. (|^)), is very interesting. There 
is another example |2^,|2^|, the off-shell massless 
ladder three- and four-point diagrams with an ar- 
bitrary number of loops, when such a generaliza- 
tion went in the CL; direction (for details, see 0). 
It could be also noted that the the two-loop non- 
planar (crossed) three-point diagram gives in this 
case the square of the one- loop function, (CI2 (#)) 2 
(cf. Eq. (23) of 0), leading to the structure 
(Cl2(7r/3)) 2 in the symmetric (pf = p 2 ) case. Re- 
cently, these constants have been also found in 
massive three- loop calculations [ p9|p6||30|l . 

The construction of analytic continuation of 
the generalized log-sine functions should be in- 
vestigated in more detail. In fact, it may require 
including some other generalizations of polyloga- 
rithms (see, e.g., in Ref. 
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